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We analyze the structure and stability of singular singly quantized vortices in a rotating spin-1 
Bose-Einstein condensate. We show that the singular vortex can be energetically stable in both 
the ferromagnetic and polar phases despite the existence of a lower-energy nonsingular coreless 
vortex in the ferromagnetic phase. The spin-1 system exhibits an energetic hierarchy of length 
scales resulting from different interaction strengths and we find that the vortex cores deform to a 
larger size determined by the characteristic length scale of the spin-dependent interaction. We show 
that in the ferromagnetic phase the resulting stable core structure, despite apparent complexity, 
can be identified as a single polar core with axially symmetric density profile which is nonvanishing 
everywhere. In the polar phase, the energetically favored core deformation leads to a splitting of 
a singly quantized vortex into a pair of half-quantum vortices that preserves the topology of the 
vortex outside the extended core region, but breaks the axial symmetry of the core. The resulting 
half-quantum vortices exhibit nonvanishing ferromagnetic cores. 
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I. INTRODUCTION 

In the textbook examples of superfiuids, liquid 4 He [l| 
and Bose-Einstein condensates (BECs) Q of spinless or 
spin-polarized atoms, quantized vortices occur as quan- 
tized circulation around an empty vortex core whose size 
is determined by a characteristic healing length. In a 
BEC of atoms whose spin degree of freedom is not frozen 
by magnetic fields Q spin rotations and condensate phase 
combine to form a larger set of physically distinguishable 
degenerate states. This is analogous to liquid 3 He where 
superfluidity is formed by Cooper pairs of fermions that 
exhibit a nonzero spin and orbital angular momentum, 
resulting in a rich phenomenology of phases with differ- 
ent broken symmetries Q. A variety of different vortex 
configurations Q and other defects and textures @ have 
been theoretically studied and experimentally observed 
in the resulting multi-component order parameter man- 
ifolds of supernuid liquid 3 He. There are obvious paral- 
lels @ to similar objects in cosmology 0] and quantum 
field theory @. 

Consequently, it is not surprising that in multicompo- 
nent BECs, there has been an increasing interest in the- 
oretical studies of topological defects and textures [9l-[39| 
as well as vector solitons (4(jl - l47l |. The theoretical work 
has been produced parallel with a rapid experimental 
progress on spinor BECs, e^g.jhi controlled preparation 
of coreless spinor vortices j48l45C| , in the studies of spin- 
texture formation (5lL HH , and in nonequilibrium vortex 
production during rapid phase transitions (53j . 

An atomic spin-1 BEC exhibits two phases of the 
ground-state manifold, ferromagnetic (FM) and polar, 
with distinct broken symmetries. In the FM phase the 
ground state of a sufficiently rapidly rotating atom cloud 
is formed by nonsingular, coreless vortices in which the 
ord er p arameter is well-defined everywhere [l?], 0, [H, 
|28|. 1.3(1 . Similar coreless vortices were first described in 



supernuid liquid 3 He [5J, [55| and were recently experi- 
mentally phase- imprinted on a spinor BEC [lU. Due to 
the topology of the FM ground-state manifold, which is 
defined by the group of spin rotations, it is also possible 
to form a singular vortex (l5l . Il6j . whose stability and 
structure, however, is much less well understood. 

Here we show that singular, singly quantized vortices 
can be energetically stable in both FM and polar phases 
of a spin-1 BEC. In the FM phase this is despite the fact 
that the coreless vortex has a lower energy. Although a 
singular vortex would also not be nucleated by rotation, 
once created, for example by phase-imprinting, it can be 
stabilized in a rotating trap. In the polar phase, the 
singular vortex undergoes a core deformation to a pair of 
half-quantum vortices in an extended vortex core region 
where the broken order parameter symmetry of the polar 
ground-state manifold is restored (see Fig.[T|). 

In a singular defect the singularity of the order param- 
eter is contained by a defect core. Unlike scalar super- 
fluids, in the spinor BEC this does not imply that the 
density must vanish: it is also possible to accommodate 
the singularity by requiring the spinor wavefunction to 
be orthogonal to the ground-state manifold at the precise 
location of the singularity. The different possibilities for 
the defect core structure leads to an energetic hierarchy 



of different characteristic length scales [23[: Depending 
on the ratio of the spin £p and density £„ healing lengths 
associated with the two interaction strengths of spin-1 
BEC, it can be energetically more favorable to force the 
order parameter value to be orthogonal to the ground- 
state manifold at the defect singularity than to force the 
density to zero. This can lead to unexpected core struc- 
tures. In Rcf. 23'j it was shown how in the polar phase 
of a spin-1 BEC a singular point defect with a vanishing 
density at the singularity can spontaneously deform to a 
ring defect with a FM core and a nonvanishing density. 
Recently symmetry classification using homotopy theory 
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FIG. 1. Schematic illustration of two vortex-core structures 
with the same topology for a singly quantized singular vortex 
in the polar phase of a spin-1 condensate. In (a) the atom 
density vanishes at the vortex-line singularity with the core 
size determined by the characteristic length scale £ n (heal- 
ing length) associated with the spin-independent interaction 
strength. In (b) the atom density is nonvanishing in the core 
region whose size is determined by the characteristic length 
scale £f of the spin-dependent interaction strength. The vor- 
tex line singularity has now split into two half-quantum vor- 
tices with the atoms in the ferromagnetic phase at the precise 
location of the singularities. In both figures we show the ne- 
matic axis as a dashed line and the dotted line in (b) indicates 
a disclination plane for the nematic axis. Inside the core re- 
gion (shaded area) of (b) the broken symmetry of the polar 
ground state manifold is restored (as explained in the text). 
Outside the core the topological properties of the vortex are 
the same as those in (a). 



was used in the analysis of defect cores in Ref . [56[ . 

In a spin-1 BEC the polar and FM phases differ by the 
local expectation value of the spin magnitude. The size 
of the filled vortex core is then determined by £p that 
defines the length scale over which the spin magnitude 
heals when locally perturbed. This is in general much 
larger than the size of a density-depleted vortex core. 
Large wavcfunction gradients close to the defect singu- 
larity result in a large order-parameter bending energy. 
Energetically, the system therefore prefers the larger core 
size and a nonvanishing atom density with correspond- 
ingly lower bending energy. Outside the filled core region 
of size £p the topology of the vortex is the same as in the 
case of a zero-density vortex line. In that region the or- 
der parameter bending energy is not sufficient to excite 
the system away from the ground-state manifold and we 
find a well-defined broken order parameter symmetry of 
either the polar or FM phase. It is only inside the filled 
core of size £p where the vortex structure differs. Inside 
the core the order parameter bending energy restores the 
order parameter symmetry of the full spin-1 condensate 
wavefunction by exciting the system out of the ground- 
state manifold by mixing the polar and FM phases. In 
our numerical simulations this is indicated by a continu- 
ously varying spin magnitude across the vortex core. 

We show that in the case of a singular singly quantized 
FM vortex the apparent complexity of the core can in this 
way be explained as the formation of a single core with 
| (F) | = at the singularity: An initial singular vortex 



is formed by overlapping vortex lines in the spinor com- 
ponents. As the system relaxes, the density depletion is 
avoided by separating the vortex lines. By a rotation of 
the spinor basis, the axial symmetry of density profiles 
of the individual spinor components with perfectly over- 
lapping vortex lines is explicitly restored, and the vortex 
is identified as one in which the spin vector winds by 
2tt around a polar core. We find both axisymmetric and 
nonaxisymmetric solutions for the singular vortex in the 
FM phase, indicating close energetic degeneracy of the 
solutions. 

In the polar phase, we use the same analysis to show 
how it is energetically favorable for the system to spon- 
taneously break axial symmetry by splitting the core of a 
singly quantized vortex into two FM cores. The resulting 
spinor wavefunction shows a complex combination of vor- 
tex lines with highly deformed anisotropic cores. How- 
ever by transforming the spinor wavcfunction to the basis 
specified by the direction of the spin in the FM cores, the 
spinor structure is identified as a pair of half-quantum 
vortices. Outside the deformed core, the topology of the 
initial singly quantized vortex is preserved. The splitting 
can be understood from the nematic symmetry proper- 
ties of the polar order parameter. A stable nonaxisym- 
metric singular vortex with a nonzero superfluid density 
at the core has been theoretically predicted [13, [HI and 
experimentally observed (59j in superfluid liquid 3 He. 

Here we analyze the energetic stability of the singu- 
lar vortices and explain the structures of the their cores 
by numerical simulations in the framework of mean-field 
theory. The paper is organized as follows: In Sec. [TT| we 
give a brief overview of mean-field theory for the spin-1 
BEC and explain the general concepts used in our anal- 
ysis. In Sec. IIIII we demonstrate the energetic stability 
and explain the core structure of a singular vortex in the 
FM phase. Section |IV] applies a similar analysis to study 
energetic stability and identify the core-structure of a sin- 
gular vortex in the polar phase. We conclude with a brief 
summary of our findings in Sec. |Vl Analytic properties 
of the vortex solutions and the basis transformations are 
provided in Appendix |A"1 

II. SPIN-1 MEAN-FIELD THEORY 

In our analysis of singular vortices of a spin-1 atomic 
BEC, we consider the classical (Gross-Pitacvskii) mean- 
field theory of a harmonically trapped system that re- 
sults in a spatially nonuniform atom density. In an op- 
tical trapping potential the atomic spin is not frozen by 
magnetic fields and the spin-1 BEC is represented by a 
normalized three-component spinor £(r) in the basis of 
spin projection onto the z axis. Together with the density 
n(r) = |*I/(r)| 2 this specifies the macroscopic condensate 
wavcfunction 

, , / C+(r) \ 

vf(r) = V^WC(r) = yfifi) Co(r) , C f C = 1. (1) 
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The Hamiltonian density in the frame rotating with fre- 
quency f2 around the z axis is [U, [60j | 



H = f\W\ 2 + V(v)n + C ^+ C ^ 
2m 2 2 



gin ( B • F ) + g 2 n ( B • F )-Q(L 



(2) 



where V(r) is an external trapping potential, m is the 
atomic mass and (L z ) = —ih^'(xd y — yd x )^ denotes the 
z-component of the angular momentum operator. The 
spin operator F, whose expectation value (F) = (^FafsCp 
yields the local spin vector, is given by a vector of spin-1 
Pauli matrices. The first two terms in the second line 
of Eq. (J2J describe linear and quadratic Zeeman shifts, 
respectively, in the presence of a weak external magnetic 
field B. In this paper we assume an axially symmetric 
harmonic confinement such that 
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V(r) = -to [uji(x 2 + y A ) + uj 1 z z 



(3) 



from which we also define the transversal oscillator length 

The spins of two colliding spin-1 atoms may combine 
to a relative angular momentum of cither or 2. This 
implies that the contact interaction results from two dif- 
ferent contributions, corresponding to the two scatter- 
ing channels with different s-wave scattering lengths ao 
and (Z2- The two scattering contributions lead to the 
two interaction terms in Eq. ([2]). The strength of the 
interactions may be calculated using angular momen- 
tum algebra and we have Co = 47rft 2 (2a2 + oo)/3to and 
c-2 = 4:7rh 2 (a2 — ao)/3m In addition the Hamiltonian 
density ^ may include magnetic dipolc-dipole interac- 
tion terms that can influence the spin textures [5ll l6l| . 

The spin-dependent interaction term c-i in Eq. @ de- 
termines the spin magnitude in a uniform ground-state 
spin distribution. If C2 < 0, as is the case for 87 Rb, 
the spin-dependent contribution to the interaction en- 
ergy will favor the FM state with |(F)| = 1 throughout 
the BEC. Conversely if C2 > 0, as for 23 Na, the polar 
state with | (F) | = will be favored. 

The two interaction strengths Cq and C2 are each as- 
sociated with a characteristic length scale. From the 
spin-independent interaction we can derive the healing 
length £ n = (87rcorj) -1 / 2 that defines the length scale 
over which the density heals around a local depletion of 
the atom density Q. This phenomenon is similar in a 
scalar BEC, which exhibits a healing length depending 
on the atom density and the scattering length. Due to 
the spin-dependent interaction term, we now, however, 
have an additional healing length, analogously given by 
£f = (87r|c2|n) -1 / 2 . This defines the length scale over 
which the spin magnitude |(F)| heals when locally per- 
turbed. 

As in a scalar BEC, single- valuedness of the order pa- 
rameter may be maintained at a defect singularity by 
requiring that the density vanishes there. The size of 



the defect core is then given by the healing length £„. 
However, the spinor order parameter makes it possible to 
maintain a nonzero density at the cost of requiring that 
the wavefunction at the singularity becomes orthogonal 
to the ground-state manifold. For example, a singularity 
in the FM manifold where |(F)| = 1 can be accommo- 
dated by having |(F)| = on the vortex line. This con- 
stitutes a local perturbation of the spin magnitude, and 
so its length scale is determined by £p, which is usually 
larger than £„. The energetic cost of the local change in 
spin magnitude due to increased interaction energy may 
be offset by the lower bending energy in the larger core. 

The two phases of the spin-1 BEC are described by 
very different order parameters, which leads to dramati- 
cally different possible vortex states. In the following we 
shall first consider the FM phase and show that a singu- 
lar, singly quantized vortex can be energetically stable, 
and describe how its core structure can be understood 
in terms of the energetics of characteristic length scales 
[Sec. IIII] . We will then apply a similar analysis to show 
how the deformed core of a stable singly quantised vor- 
tex in the polar phase can be identified as a pair of half- 
quantum vortices [Sec. IIV] , 



III. STABILITY AND CORE DEFORMATION 
OF A SINGULAR FERROMAGNETIC VORTEX 

We first consider vortices in the FM phase of a spin- 
1 BEC. The system becomes FM when the interaction 
term C2 < in the Hamiltonian @; energetically it is 
then favorable to maximize the spin magnitude every- 
where in space, so that |(F)| = 1. A general FM spinor 
wavefunction can be constructed from the representa- 
tive spinor ( = (1,0, 0) T with (F) = z by incorporat- 
ing a macroscopic condensate phase <f> and by a spin ro- 
tation U(a, P, 7) = exp(— iF z a) exp(—iF y /3) exp(— iF z j), 
defined by three Eulcr angles. We obtain 
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V2e~ ia cos 2 f 
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where 7' = 7 — 4>. The local spin vector is then given by 
(F) = (cos a sin /3, sin a sin /3, cos /3). 

Order-parameter space is the manifold of energetically 
degenerate spinors £• Degenerate FM spinors [Eq. (QJ] 
differ only by rotations in spin-space given by the Eulcr 
angles a, (3 and 7'. The order-parameter space therefore 
corresponds to the group of three-dimensional rotations 
SO(3). 

The topological stability of line defects is character- 
ized by the way closed contours encircling the defect map 
into order parameter space [62| . If the order-parameter 
space image of such a closed loop can be contracted to 
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having a corclcss vortex, which may be written 
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FIG. 2. (color online) Schematic illustrations of FM vortex 
states, (a) The non-singular, coreless vortex is formed as a 
combined disgyration of the spin vector (cones) and a spin 
rotation about the local spin vector (indicated by the orthog- 
onal vectors). The vortex is nonsingular and the spin texture 
is continuous, (b) The singular FM spin vortex is formed 
as a radially-oriented disgyration of the spin vector around 
the singular core, which is filled by the polar phase, (c) The 
class of singular vortices contain serveral different spin config- 
urations, for example the cross-disgyration shown. They can 
all be deformed to the singular spin vortex through local spin 
rotations, (d) Energies of the stable coreless and singular vor- 
tices as functions of the frequency of rotation in an isotropic 
trap using Nc = lOOOfi^x^ and Nc 2 = -320fiu±/i. The 
coreless vortex is lower in energy in the whole frequency range. 



a point, the defect is not topologically stable. SO(3) 
may be represented geometrically as S 3 (the unit sphere 
in four dimensions) with diametrically opposite points 
identified. The only closed loops that cannot be con- 
tracted to a point are those connecting such identified 
points an odd number of times (but these loops can all be 
deformed into one another) . There are therefore only two 
distinct classes of vortices: singular vortices correspond- 
ing to non-contractible loops, and non-singular vortices 
corresponding to contractible loops (j| Ho| • All the sin- 
gular vortices with an odd-integer winding number are 
therefore topologically equivalent to a singly quantized 
singular vortex and all singular vortices with an even- 
integer winding number are topologically equivalent to a 
non-singular vortex-free state. Mathematically, this is in- 
dicated by the first homotopy group of SO (3) which has 
two elements [7Ti(SO(3)) = Z2] that represent the two 
topological equivalence classes for the vortices. Typical 
examples of a non-singular coreless vortex forming a con- 
tinuous spin texture and singular spin vortex with radial 
and cross disgyrations of the spin vector are schemati- 
cally illustrated in Fig. [2] 

In the FM phase, circulation need not be quantized Q . 
A striking manifestation of this fact is the possibility of 
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where ip is the azimuthal angle and the Euler angle ft 
varies with the radial distance p = \/ x 1 + y 2 from the 
z axis. The spin texture is kept continuous by having 
f3 —> as p — >• 0. The corclcss vortex carries angu- 
lar momentum, yet can be continuously deformed into 
a vortex-free state. Because the FM order parameter 
is well-defined and non-singular everywhere in space (a 
coreless vortex forms a continuous spin texture) the core- 
less vortex belongs to the class of non-singular vortices. 
It is the spinor-BEC analogue of the Anderson- Toulouse 
and Mermin-Ho vortices in 3 He 0, [55j], which differ by 
the imposed boundary conditions at the container wall. 
The superfiuid velocity [§] 



v = (1 — cos f3)(p 

mp 



(6) 



goes smoothly to zero at the center of the vortex but 
increases away from it as j3 increases. The spin forms a 
characteristic fountain texture. In the atomic gas, there 
is no hard container wall, and the amount by which (F) 
turns as p increases from the vortex center to the edge 
of the cloud is not fixed. The total circulation can thus 
vary smoothly as the value of (3 at the edge of the cloud 
adapts to the imposed rotation [63| . 

The simplest way to construct a singly quantized sin- 
gular vortex in the FM phase is as a 27r-winding of the 
condensate phase <f>. The vortex can then be described 
by the spinor 




V2cos 2 § 
sin/3 

n2 f± 
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where the density is required to vanish on the singular 
vortex line along the z axis (where all the three spinor 
components arc singular). The Euler angle j3 is arbitrary 
but constant, giving a uniform spin distribution (which, 
without loss of generality, we assume to be in the xz 
plane such that a = in Eq. @). 

We may continuously deform £ s into another vortex in 
the same class of topological line defects through purely 
local operations. For example we may rotate the spins 
into the radial disgyration (2tt rotation) of the spin vector 
shown in Fig. HJb). This spin structure is derived from 
Eq. ((4]) by letting a — ip while 7' = 0, 
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yielding a singular vortex with a circulation of the spin 
around the core (spin vortex). This is similar to a radial 
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disgyration of the angular momentum in an analogous 
vortex structure in 3 He Q. The FM order parameter is 
still singular at p = because a singularity is introduced 
in the FM spin vector. At j3 = ir we recover the vor- 
tex in Eq. ([7]) with only one spinor component occupied. 
Further local rotations of the spin allow us to construct 
additional members of the family of singular FM vor- 
tices. If we locally rotate all spins through ir/2 around 
the z axis in Fig.[2jb) or Eq. ([8]), we change from the ra- 
dial to a tangential disgyration, where (F) = if. A spin 
vortex could also be constructed from Eq (QJ by choos- 
ing a = —ip. The radial disgyration is then replaced by 
the cross disgyration illustrated in Fig [2jc) . Because the 
SO(3) order-parameter manifold allows only two topo- 
logically distinct classes of vortices, all singly quantized, 
singular vortices can be transformed into each other by 
local spin rotations and this family of vortices is indeed 
quite large. 

In order to determine the energetic stability of the 
vortex configurations and stable vortex core structures, 
we numerically minimize the energy of specific vor- 
tex states belonging to distinct topological equivalence 
classes. The energy relaxation is done by numerically 
propagating a coupled set of Gross-Pitaevskii equations 
derived from Eq. @ in imaginary time using a split-step 
algorithm [64|. We consider an isotropic trap in a rotat- 
ing frame with the nonlinearities Aco = I000hui±lj_ and 
-640Huixlj_ < Nc 2 < -WHcj±1^, where N is the total 
number of atoms. As an initial state for a singular singly 
quantized vortex we take the vortex of Eq. ([7]) in which 
case each spinor component exhibits a singly quantized 
vortex line. These all perfectly overlap with a vanishing 
density at the core. We also perform an energy mini- 
mization of the coreless, non-singular vortex of Eq. 
The coreless vortices have been shown to exist in the 
ground state of sufficiently rapidly rotating FM spin-1 
BECs [13, El, HI HI , and increasing the rotation rate of 
a vortex-free cloud is predicted to result in nucleation of 
coreless vortices in the system. 

We find a single coreless vortex to be energetically sta- 
ble in a sufficiently rapidly rotating trap, as shown in 
the stability diagram of Fig. EJa) . Figure SJa) shows the 
characteristic fountain-like spin texture of the stable vor- 
tex. At slow rotation speeds the vortex exits the atom 
cloud and at faster rotation rates we observe nucleation of 
additional coreless vortices to the system. The threshold 
rotation frequency is increased at stronger nonlinearities. 
Our findings are consistent with those in Ref. [HI ], 

For the singular initial-state vortex the corresponding 
stability diagram is displayed in Fig. Efb). Although its 
core structure is deformed during energy relaxation (as 
we will discuss below), we find that the singular vortex 
is energetically stable for a range of rotation frequencies 
at all investigated values of c 2 - This energetic stability 
of the singular vortex seems surprising since there also 
exists a stable coreless vortex with lower energy at the 
same rotation frequencies and nonlinearities. Our nu- 
merics also show that coreless vortices will nucleate due 
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FIG. 3. (color online) Energetic stability of the coreless (a) 
and the core-deformed singular vortex (b) in an isotropic trap 
for varying spin-dependent interaction strength C2 < 0. The 
spin-independent interaction is fixed at Nco = 1000hLO±lj_. 
Blue dots (•) indicate that the vortex is energetically stable. 
A black plus (+) indicates where the initial vortex leaves the 
cloud, whereas red crosses ( x ) mark where additional vortices 
nucleate due to rotation. A blue (black) vertical line marks 
c /c 2 ~ -216 relevant for 87 Rb [|H]. Note that with the 
parameters used here, this yields 7V|c2 = The line 

thus falls very nearly on top of the vertical axis and the two 
cannot be distinguished in the figure. 
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FIG. 4. (color online) Numerically calculated spin textures 
in the stable FM vortex states in a rotating trap. The spin 
vector is shown in a cut perpendicular to the z axis (the axis of 
rotation), (a) The spin vector in the coreless vortex exhibits 
a characteristic fountain-like structure and maintains | (F) | = 
1 everywhere, (b) In the relaxed singular vortex, the spin 
vector winds by 2-n around the x axis on a path encircling the 
singular vortex core (indicated by the dot), in which |(F)| — > 
0. This texture can be continuously deformed into that shown 
in Fig. Ob). 



to rotation, whereas singular vortices will not. A com- 
parison between the numerically calculated energies of a 
stable coreless vortex and a stable singular vortex as a 
function of the rotation frequency in an isotropic trap and 
with Ac = lOOOfa^l, Nc 2 = -320fto;_L^i is shown in 
Fig. [2Id). This may be contrasted with the vortex ener- 
getics of coreless and singular vortices in supcrfiuid liquid 
3 He-A, where the singular vortex has lower energy, but 
the energy barrier for nucleation of the singular core is 
higher than that for forming a coreless vortex [66| . Singu- 
lar vortices can be created by cooling a rotating normal 
fluid through the supcrfiuid transition. 

The coreless and the singular vortices belong to dis- 
tinct topological equivalence classes and they cannot be 
continuously deformed to each other. For the singular 
vortex to decay, the rotation frequency has to be suffi- 
ciently slow so that the vortex can exit the atom cloud 



and be replaced by a nucleating coreless vortex that en- 
ters from the edge of the cloud. We find a range of 
frequencies and nonlincarities [Fig. EJb)] for which the 
singular vortex remains in the atom cloud and no addi- 
tional coreless vortices nucleate. A single, singly quan- 
tized singular vortex thus represents a local minimum of 
the energy, topologically protected against decay to the 
lower-energy coreless vortex. 

After demonstrating that the singular singly quantized 
vortex of Eq. ([7]) is energetically stable, we next study 
its vortex core structure after the energy relaxation. The 
resulting vortex configuration with a stable vortex core 
is shown in Fig. [5ja) . The vortex lines in the different 
spinor wavefunction components have moved apart and 
no longer spatially overlap. We show in Fig. EJa) a ID 
density cut along which the spatially separated vortices 
are aligned. The vortex line of the Co component is lo- 
cated at the center of the trap and the vortices of the £± 
components are symmetrically displaced from the center. 
This split-core solution appears to break the explicit axial 
symmetry of the spinor component densities in Eq. ([7]). 
A similar core splitting has previously been demonstrated 
in 2D numerical simulations in Ref. (HI- We will show 
below how it is beneficial to analyze the vortex core us- 
ing a spinor basis transformation. In particular, after 
an appropriate transformation we can easily identify the 
location of the vortex, the nonvanishing atom density 
at the vortex line singularity, and the axially symmetric 
density profiles of the spinor components in the new ba- 
sis representation. In the vortex configuration displayed 
in Fig. E][a) we may then identify the split-core vortex as 
spin winding around a core of nonvanishing atom density. 

In order to analyze the vortex configuration of Fig.[5ja) 
we perform a basis transformation for the spinor wave- 
function. We transform the split-core spinor to the ba- 
sis where spin is quantized along the x axis as C^ x > = 
C/ _1 (0, 7r/2, 0)C^\ explicitly indicating the spinor basis 
by superscripts. In £W the vortex appears as an oppo- 

(x) 

site winding of the phase in the two components Q_ . 
These vortex lines again overlap as shown in Figs. Efb) 

and[B{b). Crucially, there is no vortex line in (q , an( i 
this component therefore fills the vortex cores of the two 
other components so that the density is nonvanishing ev- 
erywhere. The single vortex core, which is readily appar- 
ent from Fig. [51 is thus explicitly restored in Q x > by the 
transformation to the 'natural basis' of the vortex. We 
identify the spinor wavefunction resulting from the basis 
transformation now as having the same structure as the 
singular vortex, defined in Eq. . In Appendix [A] we 
show how the core deformation and the relation between 
and C (2) can be understood qualitatively through an 
analytic treatment. 

The spin structure of the stable vortex is shown in 
Fig. 0|b). The vortex line is oriented along the z 
direction — the axis about which the trap is rotating. 
However, the spinor takes the form of Eq. © in the ba- 
sis defined along the (co-rotating) x axis. This vortex 



6 




FIG. 5. (color online) Split core of the singular FM vortex 
and restoration of a single core with explicit axial symmetry 
of the spinor-component densities illustrated by isosurfaces 
of the spinor wavefunction components n^+| 2 (red/medium 
gray), n |Co| (green/light gray) and n |£_| (blue/dark gray), 
(a) In the spinor basis along the z axis, the vortex lines in 
C+> Co Z ' an d C-* 1 separate and the atom density is nonzero 
everywhere, (b) The axial symmetry of the density in each 
spinor component is restored by transforming the spinor to 
the basis of spin projection onto the x axis. Vortex lines with 
opposite circulation in overlap. Q B ' (not shown) does not 
exhibit any vortex line (cf. Fig- Ob)). See also Appendix [Al 
for a qualitative analytic discussion of the relation between 
C w and ( {z) . 




FIG. 6. (color online) (a) Densities in the three spinor com- 
ponents (red line marked by +), Q (green line marked 
by 0) and (blue line marked by — ) on the axis connect- 
ing the vortex lines in the spinor components (cf. Fig(5ja)). 
(b) Densities in Qf^ (red line marked by +), Q x ' (green line 
marked by 0) and f_ (blue line marked by — ) on the same 
spatial axis after spinor basis transformation. |(F)| (black 
dash-dotted line) goes to zero in the vortex core (the apparent 
nonzero minimum is due to finite numerical resolution) which 

(x) 

is filled by Q i keeping the density nonzero everywhere. 

is singular, preserving the topology, and can be reached 
from Eq. {7} by local spin rotations and could similarly 
be continuously transformed into the singular spin vor- 
tex [Fig. EJb)]. The stable vortex core Fig. |4|b) has a 
broken spatial parity (the spin profile has an antisym- 
metric spatial parity close to the vortex core). This spin 
profile is nonaxisymmetric. We also find a stable axisym- 
metric vortex core. This is achieved by starting energy 
relaxation from Eq. ([SJ, such that the radial disgyration 
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FIG. 7. (color online) Axially symmetric spin vortex, (a) 
Densities of the spinor components together with the spin 
magnitude along a radial cut (lines and labels as in Fig. EJa)). 
The vortex lines in £± overlap perfectly at the position of the 
vortex core, (b) Relaxed spin profile in the xy plane, showing 
the characteristic radial disgyration of the spin vector around 
the singular core. At large radii the spin vector bends out of 
the xy plane. The vortex line singularity is marked by a dot 
at the center. 




FIG. 8. (color online) Spin magnitude |(F)|, showing the core 
of the stable singular vortex of Fig.[S]in the FM phase, shown 
in the xy plane (a) and yz plane (b). Outside the vortex core, 
|(F) | = 1 (dark red/dark gray) in the FM order parameter 
manifold. In the core, the singularity is accommodated by 
enforcing | (F) | = (white) while maintaining nonzero density. 
The size of the core region is determined by the spin healing 
length £f- 



of the spin vector is present already in the initial state. 
The spinor components and the resulting spin profile are 
shown in Fig. [JJ The dependence of the final configura- 
tion on the initial state indicates a close energetic degen- 
eracy of the two solutions. 

To understand the vortex core deformation it is also 
beneficial to compare the initial-state singular vortex of 
Eq. ([7]) to the vortex obtained in the energy minimiza- 
tion. In Eq. (J7J each spinor component exhibits a singly 
quantized vortex. These overlapping vortex lines imply 
that the total density n(r) must be zero on the singular 
line in order to maintain single-valuedness of the order 
parameter. The size of the vortex core is determined 
by the healing length £„. The density depletion can be 
avoided by splitting the vortex core such that the vor- 
tex lines in the spinor components no longer overlap. 



Since the total condensate density then does not vanish 
at the vortex line where the order parameter is singu- 
lar, we must now require that the spinor wavcf unction 
becomes orthogonal to the ground-state manifold at the 
vortex singularity. In the FM manifold | (F) | = 1, so at 
the vortex line we must have | (F) | = 0, which represents 
the spin magnitude of the polar phase. 

The spin magnitude of the numerically calculated sin- 
gular vortex core is displayed in Figs. Eljb) and |U We 
find that the value of the spin magnitude indeed rapidly 
approaches zero close to the vortex line singularity (the 
small deviations from zero are due to spatial resolution 
of the numerics). This indicates the formation of a po- 
lar vortex core, constituting a local violation of the spin 
condition for the ground-state FM manifold. An ana- 
lytic description of the vortex solution is provided in Ap- 
pendix [XJ The size of the vortex core is determined by 
the spin healing length £p. The splitting is then ener- 
getically favorable when £p allows a larger core size (i.e., 
when £f > £„) such that the energy cost of violating 
| (F) | = 1 is smaller than that of depicting the density. 

We find that the region where the spin magnitude de- 
viates from |(F) | = 1 extends over the entire core size, 
determined by the spin healing length £p. Outside the 
core region of the vortex the symmetry of the spin-1 
BEC is broken according to the FM energy condition of 
the spin-dependent interaction energy, so that we have 
|(F) | = 1. Close to the singular vortex, however, the 
order-parameter bending energy restores the symmetry 
of the full spin-1 BEC wavefunction (S 5 determined by 
a normalized spinor wavefunction of three complex com- 
ponents), mixing the FM and polar phases. The bending 
energy is enhanced very close to the vortex singularity 
due to the large density gradient contributions that ex- 
cite the system from the FM ground-state manifold. An 
analogous core deformation was previously found for a 
singular point defect in a polar spin-1 BEC in Ref. [23| . 
In that case an isotropic point defect with a vanishing 
density deformed to a ring defect with a FM core. This 
effect is closely related to the deformation of the core of 
a singular vortex in the polar phase described in Sec. IIVI 

In experiments a stable singular vortex could be pre- 
pared in a controlled way by phase-imprinting the initial 
singular vortex state of Eq. (J7J in a rotating trap, so that 
the parameter values of the system belong to the stable 
region of the stability diagram displayed in Fig. [3f b) . 
The initial-state vortex [Eq. (J7J)] is composed by per- 
fectly overlapping singly quantized vortices in each of the 
spinor components. These could be phase-imprinted us- 
ing previously realized experimental techniques (67rlo9| . 
The stability diagram also indicates the conditions under 
which a singular vortex created in a phase transition (53j 
could potentially be stabilized. 

In the above analysis, we have allowed the magnetiza- 
tion M z = N + — , where N± are the total populations 
of £± , to vary during the relaxation process. This in prin- 
ciple allows a spontaneous magnetization to develop in 
the system. In experiments, dissipativc relaxation of en- 
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FIG. 9. (color online) Stability of the FM coreless (a) and 
singular (b) vortices in a highly oblate trap, u)±/ui z — 0.1, 
with Nco = 50hw±l 3 ± and -4&f\w±l\ < Nc 2 < -5?iu±l±. 
(Symbols as in Fig. [3]) 



ergy due to atomic collisions may frequently conserve the 
magnetization. We have therefore also performed calcu- 
lations where a weak magnetization is enforced through- 
out the energy-minimization procedure. We find the re- 
laxed vortex configurations of Figs. [5][H however, quali- 
tatively similar to the ones where no fixed magnetization 
was forced. 

Thus far we have considered an isotropic trap. We 
find that the results are qualitatively similar in an oblate 
trap with uj±/lu z =0.1. We find that also in this regime, 
the singular vortex represents a local energetic minimum 
and is stable for a range of f2, again despite the fact that 
a lower-energy coreless-vortex solution exists. The pa- 
rameter regions allowing stable nonsingular coreless and 
singular split-core vortices in the oblate trap are shown in 
Fig. [9] In the oblate trap case the simulations with fixed 
magnetization produce more symmetric vortex configu- 
rations than in the simulations where the magnetization 
varies freely but the qualitative features are the same in 
the two cases. 

Applying a weak external magnetic field introduces 
a Zeeman shift between the spinor components accord- 
ing to Eq. ([2]). The spinor nature of the BEC is re- 
tained as long as the applied field is not too strong 
<7i |B | , <72|B| 2 < fi, where fi denotes the chemical poten- 
tial. In the case of a small linear Zeeman splitting <?i|B| 
(taking B along the z axis) we find that the coreless and 
singular vortices are both stable, with the coreless vor- 
tex lower in energy. The Zeeman splitting will tend to 
align the spins with the applied field. This causes the 
energy of the coreless vortex to increase as maintaining 
the fountain-like spin structure becomes energetically less 
favorable. Thus we find that for <?i|B| > Q.2huj± the core- 
less vortex is no longer stable. The singular vortex, on 
the other hand, remains energetically stable for all <7i|B| 
considered (up to 0.8hu±), as shown in Figs. ITUT a) and 
(b). For a sufficiently large linear Zeeman splitting the 
ideal spinor basis to analyze the singular vortex core be- 
comes the one defined by the magnetic field. 

A quadratic Zeeman splitting, on the other hand, does 
not destroy the stability of the coreless vortex, but for 
<72|B| 2 < — 0.1hu± the singular vortex is no longer en- 
ergetically stable [Fig. [TUtc) and (d)]. For a sufficiently 
large positive quadratic Zeeman splitting the ideal spinor 
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FIG. 10. (color online) Effects of linear and quadratic Zee- 
man splitting on the stability of the FM vortices in an oblate 
trap (uj±/lj~ = 0.1). In all panels Nco = 50faj±Zx an d 
Nc2 — — 10/ojx'x- ( a ) Stability of a coreless vortex in the 
presence of linear Zeeman splitting. The vortex state be- 
comes unstable for gi |B| > 0.2ftu)±. (b) The singular vortex 
remains stable despite linear Zeeman splitting, (c) Stability 
of the coreless vortex in the presence of quadratic Zeeman 
splitting. A stable region is found at all investigated values 
of (72 1 33 1 . (d) Stability of the singular vortex in the presence 
of quadratic Zeeman splitting. The vortex becomes unstable 
for a relatively small negative g 2 |B|. (Symbols as in Fig. O) 



basis to analyze the singular vortex core is oriented per- 
pendicular to the magnetic field. 



IV. STABILITY AND CORE STRUCTURE OF A 
POLAR VORTEX 



We now turn our attention to the polar phase of a 
spin-1 BEC. In this case the interaction term C2 > in 
the Hamiltonian @, and it is energetically favorable to 
minimize the spin magnitude everywhere in space so that 
|(F) | = 0. We take a representative spinor ( = (0, 1, 0) T 
whose macroscopic condensate spin quantization axis is 
oriented along the z axis. The general spinor wavefunc- 
tion may then be constructed from the macroscopic con- 
densate phase (j) an d the spin rotations defined by the 
Euler angles (a, f3, 7) as 





e»U(a,p,i) ( 1 





-e M sin/3 
VZcos/3 
e la sin /? 



(9) 



It is beneficial to introduce the unit vector d = 
(cos a sin /?, sin a sin (3, cos /?) that defines the local direc- 
tion of the condensate spin quantization. We may then 
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write the spinor wavefvmction in terms of d as [22 



e i<f> 

71 



—d x + idy 
V2d z 
+ idy 



(10) 



The unit vector d takes values on a sphere and the con- 
densate phase (j> on a unit circle. The state of the spinor 
wavefunction, however, remains unchanged when a n- 
rotation of (f> is combined with inverting the d vector, 
so that the states C p (</ ) ? d) = C p (0 + ""j - d) are iden- 
tical. These states must be identified to avoid dou- 
ble counting, and the order-parameter space is therefore 
(U(l) x S 2 )/1j2, from the condensate phase and rotations 
of d, factorized by the discrete two-element group Z2 due 
to the identification. The vector d is thus taken to be un- 
oriented and defines a nematic axis [2~i| . 

A singly quantized vortex in the polar phase can be 
formed as a 27r-winding of the condensate phase <f> around 
a closed loop encircling the vortex core. Choosing the 
vortex line along the z axis, we obtain 



tlp . -e ta sm/3 
C 1 = -= a/2cos/? 

e la sin /3 



(11) 



However, the nematic order also allows the formation of 
a vortex carrying half a quantum of circulation (llj (in 
an analogy to half-quantum vortices in superfluid liquid 
3 He (H), constructed as a n- winding of the macroscopic 
condensate phase together with ad-> d rotation of the 
nematic axis around a closed loop encircling the vortex 
core. For example, we may have 



-1/2 _ 




V2 



(12) 



*v 



Circulation is thus quantized in units of n, half the cir- 
culation of Eq. (fnj) . This will be a crucial observation 
when we now analyze possible deformations of the core 
of a singly quantized vortex as energy is minimized. 

In order to investigate the energetic stability of a singly 
quantized singular vortex in the polar phase of a spin-1 
BEC we numerically minimize the energy of the system 
in a rotating frame. We follow the same procedure as in 
the FM case and this time take a singular polar vortex 
of Eq. (|lip with f3 — 7r/4 and a = as an initial state of 
the numerical relaxation. Similarly as in the FM case of 
Eq (|7|), the initial state is formed by overlapping vortex 
lines in all the three spinor components. Upon minimiz- 
ing the energy, the vortex cores of the individual spinor 
components separate. However, compared with the FM 
case, the splitting is now more complicated, as shown in 
Fig. [TT] The result is highly deformed anisotropic vor- 
tex cores in the spinor components. The vortices in C+ 
and C- overlap, but the one in Co is displaced from the 
other two. There are no simultaneous density minima in 
all three spinor components, and the density is therefore 
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FIG. 11. (color online) Stable core structure of the singular 
vortex in the polar phase shown in the xy plane, (a) and 
(b): Densities in Qf' and C} "\ respectively, (c) and (d): The 

corresponding phases. (Ci z ' is identical to up to a global 
7r phase shift.) The spinor wavefunction exhibits vortex lines 
with highly deformed anisotropic cores in the spinor compo- 
nents. 



nonzero everywhere. Similar split-core solutions found 
by numerical calculation in a rotating 2D system (Io| 
have resulted in some controversy regarding the number 
vortices in the individual spinor components in the fi- 
nal configuration [l6|. In the previous 2D studies the 
stable core structures were not classified. Here we show 
now how the split core in Fig. [TT] can be identified as 
a topology-preserving splitting of the singly quantized 
vortex into a pair of half-quantum vortices as illustrated 
schematically in Fig. [TJ 

In the numerical simulations the initial state of a singly 
quantized singular vortex in Eq. ([lip is composed of three 
perfectly overlapping vortex lines in each of the three 
spinor components. The polar vortex consequently has 
a vanishing density at the line singularity of the polar 
order parameter of the spin-1 BEC. The singular vor- 
tex with zero density is energetically unstable with re- 
spect to core deformation. As the vortices of the indi- 
vidual spinor components move apart during energy re- 
laxation, the density becomes nonvanishing everywhere 
in the vortex-core region. Similarly to the FM vortex 
case, we must therefore require that the spinor wavefunc- 
tion becomes orthogonal to the ground-state manifold at 
the vortex singularity. This indicates that we must have 
I (F) I = 1 on the vortex line. We show in Fig. [T2] the 
numerically calculated vortex core structure of a stable 
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FIG. 12. (color online) Splitting of the singly quantized 
vortex into two half-quantum vortices, (a) Spin magnitude 
|(F) | [color map from white (|(F)| = 0) to red/dark gray 
(|(F) | = 1)] together with the spin vector (arrows), showing 
the FM cores with nonvanishing density. The spins are an- 
tiparallel in the two cores, (b) Nematic axis d(r) together with 
the vortex cores (indicated by green/light gray isosurfaces of 
|(F) |, with increasing spin magnitude indicated by the color 
gradient inside). Away from the vortex cores the topology of 
the initial singly quantized vortex is preserved. In the core 
region, d winds by n about each half-quantum vortex core. 
For visualization purposes, the unoriented d-field is shown as 
cones. Here a quadratic Zeeman shift has been introduced to 
ensure that d lies in the xy plane and the spins align with the 
z axis. 




vortex whose initial state is the singular singly quantized 
vortex of Eq. (jlip . The displayed spin magnitude ex- 
hibits two clearly separated cores in which the peak value 
increases to |(F)| = 1, indicating the emergence of a FM 
core region for the vortex. 

The formation of the FM cores can be understood from 
the same argument used to understand the polar core of 
the singular FM vortex in Sec. IIIII and is illustrated in 
Fig. [TJ The singular polar vortex (fTTj) , which is used as 
an initial state in the energy relaxation, implies a density- 
depleted core whose size is determined by £„. However, 
accommodating a singularity of the polar order parame- 
ter by requiring | (F) | = 1 at the vortex line means that 
the length scale, and thus the associated bending energy, 
is determined by £f- The energy of Eq. (fTTj) can thus be 
lowered by having a nonvanishing atom density and by 
extending the core size from £„ to £p. In the case of a 
polar vortex this is achieved by spontaneously breaking 
the axial symmetry and forming two FM cores by the 
mechanism sketched in Fig. [TJ The separation between 
the cores is of the order of £p, depending also on the 
angular momentum of the system when it adjusts to the 
rotation frequency and on the density gradient due to the 
trap. 

We may analyze this symmetry breaking of the vor- 
tex core by means of a basis transformation. We write 
the spinor in the basis of spin projection onto the axis 
given by the spin vector in the FM core. For the 
case of Fig. [TTTa'). the spins in the two cores align 
(antialign) with the y axis and we calculate Q v ' = 
C/- 1 ( 7 r/2,7r/2,0)C (z) . The resulting spinor then shows 
displaced vortex lines in Q?^ while the density vanishes 
in [Fig.[T3]. In the new spinor basis, the vortex lines 



FIG. 13. (color online) Spinor wavefunction of the stable sin- 
gular vortex state from Fig. [TTJ after spinor basis transforma- 
tion such that spin is quantized along the y axis, (a) and 
(b): Densities in Q?' and respectively, (c) and (d): The 

corresponding phases. The component (not shown) is 
unpopulated. The previously complex structure can now be 
identified as a pair of half-quantum vortices. 
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FIG. 14. (color online) (a) Density profiles of the spinor 
components on the axis connecting their density minima (cf. 
Fig. 1111) . Lines labeled with spinor component index. Note 
that |£±| exactly overlap, (b) Spinor-component density pro- 
files in after basis transformation (cf. Fig. Ill[) plotted 
along the axis connecting the half-quantum vortices. The 
spin magnitude |(F)| (black dash-dotted line) shows the FM 
cores. Unpopulated is not shown. 



in coincide precisely with the spin maxima, as shown 
in Fig. EH 

We can now identify the core structure emerging from 
the splitting of the singular vortex by comparing the spin- 
rotated state C^* 1 with Eq. (fT2|) . We then find that each 
vortex line in has exactly the form of a half-quantum 
vortex. The split-core configuration may thus be inter- 
preted as a splitting of the singly quantized vortex into a 
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pair of half-quantum vortices with FM cores. The topo- 
logical charges of vortices are here additive and topology 
is therefore preserved when the singly quantized vortex 
splits into the pair of half-quantum vortices. This can 
also be inferred from the behavior of the nematic axis 
d. Figure [T2jb) shows d in a numerical solution together 
with the FM cores of the half-quantum vortices. Away 
from the vortices, there is no net winding in d on a path 
enclosing the vortices. However on a path that encir- 
cles only one vortex core, d turns by n, indicating the 
emergence of a disclination plane as indicated in Fig. [T] 

As in the case of a FM vortex, the core deformation 
can be explained in terms of the vortex topology and the 
energetic hierarchy of different length scales (see Fig. [TJ . 
Outside the vortex core region of size £p, where the order 
parameter bending energy is not sufficient to excite the 
system away from the polar ground-state manifold, we 
have | (F) | = and the topological properties of the ini- 
tial singly quantized singular vortex are preserved. This 
is indicated by the unit winding of the macroscopic con- 
densate phase around any closed loop encircling the en- 
tire vortex core and by the nematic vector field outside 
the core region. It is only inside the core of size £p that 
the strong order parameter bending energy restores the 
symmetry of the full spin-1 condensate wave function by 
exciting the system out of the polar ground-state man- 
ifold and by allowing the complete range of spin values 
|(F) | from to 1. The local deformation of the core is 
topologically possible due to the nematic order of the 
polar phase, where the axis d is unoriented, with the 
opposite orientations d = — d identified. The core de- 
formation mechanism of the vortex line is related to the 
point defect deformation into a singular ring where the 
nematic order allows the spontaneous breaking of the 
spherical defect core symmetry (23|. In the £>-phase of 
supcrfluid liquid 3 He, a stable nonaxisymmetric singular 
vortex with a nonzero superfiuid density at the core was 
theoretically predicted in Refs. [13, HH and experimen- 
tally observed in Ref. [59|. The 3 He A- phase core was 
explained to consist of two half-quantum vortices. In the 
high-pressure regime the axial symmetry of the vortex is 
restored but the core can still remain in the A-phase with 
a nonvanishing supcrfluid density [70j . 

We find that the singular vortex splits into a pair of 
half-quantum vortices by the mechanism described above 
for all investigated parameter regimes. However, we find 
a critical rotation frequency of fi ~ 0.3uj± below which 
the vortices start exiting the atom cloud. FigurefTSlshows 
the energetic stability of the half-quantum vortex pair ob- 
tained from splitting of a singly quantized vortex in both 
isotropic and oblate (uj±/lu z = 0.1) traps. The energetic 
ground state of a rotating polar spin-1 BEC can consist 
of half-quantum vortices [26|, [28[, and increasing the rota- 
tion frequency leads to nucleation of more half-quantum 
vortices in addition to the split core of the initial singular 
vortex. 

The splitting mechanism of the singly quantized vor- 
tex is qualitatively similar also when a weak Zeeman 
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FIG. 15. (color online) Energetic stability of the split-core sin- 
gular vortex in the polar phase, (a) Stability in the isotropic 
trap for varying C2 using Ncq = 1000/iWjJj. The vertical 
line marks the value c /c 2 ~ 31 for 23 Na [7l[. (b) Stabil- 
ity in the oblate trap (wi/w z = 0.1) for varying ci using 
TV co = 50&Jx^x- (c) and (d): Stability of the split-core sin- 
gular vortex in the presence of linear and quadratic Zeeman 
splitting, respectively. Increasing linear Zeeman splitting ren- 
ders the singular vortex unstable, whereas the stability is ro- 
bust against quadratic Zeeman splitting. In the slowly ro- 
tating region for all panels, the instability of the split-core 
singular vortex may be either towards the vortex-free state, 
or towards the state with a single half-quantum vortex. (Sym- 
bols as in Fig. [3]) 



splitting due to a magnetic field is introduced. How- 
ever, as shown in Fig. I15f ch a linear Zeeman splitting of 
g\B > 0Ahu± causes the resulting pair of half-quantum 
vortices to become energetically unstable at all rotational 
frequencies. By contrast, the vortex pair remains stable 
above f2 ~ 0.3ujj_ for the entire range of quadratic Zee- 
man splittings considered (— 0.8fruj± < g2\B\ 2 < 0.8huj±) 

[Fig.rnd)]. 



V. CONCLUSIONS 

We have demonstrated that a singular, singly quan- 
tized vortex can be energetically stable in the FM phase 
of a spin-1 BEC, despite the existence of a stable core- 
less vortex with lower energy (increasing the rotation fre- 
quency in a vortex-free BEC also leads to nucleation of 
coreless vortices). This implies that even though singu- 
lar vortices would not be nucleated by rotation alone, a 
singly quantized vortex created, for example, by phase- 
imprinting would remain stable in the rotating system. 
This provides an interesting opportunity for controlled 
studies of a singular vortex line in a ground-state man- 
ifold with a broken SO (3) symmetry. Such a system 
supports only two topological classes of vortices: those 
that can be locally deformed to a vortex-free configu- 
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ration and those that are topologically equivalent to a 
singly quantized singular vortex. Experimentally, one 
could phase-imprint overlapping vortex lines in each of 
the three spinor wavefunction components. The result- 
ing structure represents a singular spin-1 vortex with a 
vanishing density at the vortex line. The core of such a 
vortex then deforms to a energy-minimized configuration 
within the same topological equivalence class. 

The stable vortex core in the FM phase is formed by 
nonoverlapping vortex lines in the three spinor compo- 
nents. We have demonstrated that this seemingly com- 
plex core structure can be understood in terms of the 
combination of the vortex topology and the energetics 
of characteristic length scales. By deforming the core of 
a singly quantized, singular vortex in the FM phase so 
as to maintain a nonzero density everywhere, instead ac- 
commodating the singularity by forcing |(F)| = 0, the 
gradient contribution to the energy is lowered. The rea- 
son is that the size of the defect core is then determined 
by the spin healing length £p which is in general larger 
than the characteristic size of a defect core where the 
density goes to zero. In other words, in the larger core 
size case with a nonvanishing atom density, the gradient 
energy restores the full symmetry of the spin-1 conden- 
sate wavefunction within the core region. The system 
then simultaneously exhibits two different order parame- 
ter symmetries: maximal unbroken symmetry inside the 
core of size £f and a broken symmetry (of the FM phase) 
outside the vortex core. 

The core deformation mechanism results in a singular 
vortex whose core is also filled with atoms in the polar 
phase. The spin vector winds by 2tt as the core is encir- 
cled. The single vortex core can be explicitly restored in 
the spinor by judicious choice of spinor basis. 

In the polar phase, we have shown that a singly quan- 
tized vortex is stabilized by a spontaneous breaking of 
axial symmetry. The resulting stable defect is a pair of 
half-quantum vortices with FM cores, which is stable in 
a sufficiently rapidly rotating trap. The formation of the 
FM cores avoids depleting the density in the vortex core. 
This is energetically favorable by the same reasoning that 
was applied to explain the polar core of the singular vor- 
tex in the FM phase. The resulting spinor wavefunction 
is analyzed and the vortex structure identified through a 
rotation of the spinor basis, so that in the rotated basis 
the half-quantum vortices appear as separate vortex lines 
in the £± components. 
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Appendix A: Basis transformation for FM vortex 

In the numerical simulations we found that the singu- 
lar FM vortex relaxes to a stable configuration formed 
by nonoverlapping vortex lines in the three spinor com- 
ponents, as shown in Fig.[5ja). By an appropriately cho- 
sen basis transformation we showed that this seemingly 
complex vortex structure can be identified as a single, 
singular vortex with the line singularity populated by 
atoms in the polar phase [Fig.[5jb)]. In this Appendix we 
demonstrate this basis transformation through a qualita- 
tive analytic treatment and show how the core structure 
of the singular FM vortex may be identified. 

For simplicity, we implement the basis transforma- 
tion by starting from the final configuration of the sin- 
gular vortex with a single vortex line as in Figs. [SJb) 
and [7] and rotating to the configuration of nonoverlap- 
ping vortex lines in the three spinor components [see, 
e.g., Fig. [5][a)]. The analytic expressions become notably 
simpler in the case of the axisymmetric vortex of Fig. [7J 
than with the one displaying a more complex spin rota- 
tion in Figs. ??(b) and [5j but the basic principle of the 
transformation is the same in both cases. In order to de- 
scribe the vortex of Fig. [JJwe rewrite the singular vortex 
displaying a radial disgyration of the spin vector £ sv of 
Eq. in the following form 

-, / \/2(cos</3 — i sin (/j)/ cos 2 f \ 
C=-= gsm(3 . (Al) 

^ l \ ^(cos^ + isin^/sin 2 | / 

Here we have introduced the profile function f(x,y) and 
the notation 

g sin /3 = y/2 - f 2 (1 + cos 2 0) . (A2) 

The profile function / describes the mixing of the FM 
and the polar phases in the core region of the vortex. For 
(5 7^ 7r we obtain a radial disgyration with nonvanishing 
density at the singularity (for ft = n and / = 1 we recover 
a vortex with zero density at the singularity) . We assume 
/ to be monotonically increasing from at the z axis (the 
vortex singularity) reaching 1 outside the vortex core of 
size £p . For the numerically minimized stable solution 
of Fig. [7J the parameter (3 is not constant; close to the 
vortex core we have j3 ~ 7r/2 and far away from the 
vortex /3 — > n. The magnitude of the spin vector may be 
evaluated from Eq. ([AT]) , yielding |(F)| = fy/2 - f 2 . At 
the line singularity | (F) | = representing the polar phase 
and outside the core region |(F)| = 1 that corresponds 
to the FM phase. In between |(F)| continuously varies 
between these two values, indicating the mixing of the 
two phases. 

When we perform the rotation of Eq. (|Alj) by the angle 
of — 7r/2 with respect to the x axis, we obtain the spinor 
wavefunction 

Y I I ( cos f ~~ * sm f cos P) + 9 sm P \ 
C { ~ x) = - v^/(-cos</?cos£ + £sin^) . (A3) 
\ / (cos Lp — i sin ip cos (3) — g sin /3 J 
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The spinor wavefunctions are of the form (x— Xo)+if]{y— 
Do), indicating that a singly quantized vortex line is lo- 
cated at (xq,i/q). The anisotropy of the vortex core is 
described by the parameter ?/. The singularity in Q 



therefore is on the z axis, while those in Q_ are dis- 
placed to (t^OiZ/o = 0). Here Xq is determined as 
the point at which /(xo,0) = g(x 0l 0) sin/3, resulting in 
/(xo,0) = \/2/(2 + cos 2 (3). The vortex configuration of 
Eq. (|A3|) with three spatially separated vortex lines is 
analogous to that shown in Fig. Eta). 
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